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P_i ' If expanding and contracting regions coexist in the universe, the speed of the cosmic vol- 

ume expansion can be accelerated. We construct simple inhomogeneous dust-filled universe 
models in which the speed of the cosmic volume expansion is accelerated for finite periods. 
These models are constructed by removing spherical domains from the Einstein-de Sitter 
universe and filling each domain with a Lemaitre-Tolman-Bondi dust sphere possessing the 
same gravitational mass as the removed region. This represents an exact solution of the 
Einstein equations. We find that acceleration of the cosmic volume expansion is realized in 
some cases when the size of the contracting region is comparable to the horizon radius of the 
Einstein-de Sitter universe though this model is very different from the universe observed 
today. This result implies that non-linear general relativistic effects of inhomogeneities are 
very important to realize the acceleration of the cosmic volume expansion. 
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§1. Introduction 
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The accelerating expansion of the universe indicated by observational data for 
the luminosity distance of Type la supernovae and the cosmic microwave background 
radiation is a great mystery in modern cosmology. The acceleration of a homogeneous 
and isotropic Friedmann-Robertson- Walker (FRW) universe implies the existence 
of a dark energy component of matter fields which violates the energy conditions 
and whose nature is unknown. Conversely, if the universe is not homogeneous and 
isotropic, the observational data do not necessarily indicate an accelerating expansion 
of the universe, or even if the cosmic expansion is accelerating, it does not necessarily 
imply the existence of dark energy. Thus, to account for the observational data 
without introducing dark energy, arguments concerning the effects of inhomogeneities 
have been made. Roughly speaking, there are two such arguments. One is that 
the apparent acceleration of our universe can be regarded as a result of an almost 
spherically symmetric but inhomogeneous peculiar velocity field, assuming that we 
are located in the vicinity of its symmetry center. 1 ) -4 ) With this argument, the 
acceleration of the cosmic volume expansion is not necessary. The other argument is 
that the apparent acceleration of our universe results from backreaction effects due 
to inhomogeneities in the background FRW universe. 5 )' 6 - 1 In spite of its viability, the 
former argument is simple and clear. By contrast, at present it is unclear whether 
the backreaction effects of inhomogeneities can actually accelerate the cosmic volume 
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expansion and, further, account for the observed distance-redshift relation. 7 ^ 9 ^ 

Perturbative analysis of backreaction effects on the background FRW universe 
shows that spatial inhomogeneities behave as a positive spatial curvature, and their 
presence reduces the expansion rate of the FRW universe. 10 )~ 12 ) Further, recently 
Kasai et al. conjectured that a non-linear backreaction cannot accelerate the cosmic 
volume expansion even if the inhomogeneities are very large. 13 ' By contrast, one of 
the present authors, Nambu, and a collaborator, Tanimoto, pointed out the possi- 
bility that the non-perturbative features of inhomogeneities are necessary to realize 
an accelerating expansion of the universe through backreaction effects. 14 ) They pro- 
posed a model of the universe containing both expanding and contracting regions 
and showed that a spatially averaged scale factor defined in terms of the volume of a 
spatial domain can exhibit accelerated expansion if the size of each region is properly 
chosen. 

At present, it is still unclear whether the acceleration of the cosmic volume 
expansion could be realized through the effects of inhomogeneities in the universe 
without introducing dark energy. Thus this point must be investigated before we 
attempt to explain the observed distance-redshift relation with this approach. As 
Nambu and Tanimoto have not presented a correct example of an accelerating uni- 
verse, it is necessary to determine whether their mechanism actually works. In this 
paper, in order to do this, we study the comoving volume of a highly inhomogeneous 
dust-filled universe. 

This paper is organized as follows. In §2, we give a qualitative description of 
a scenario in which inhomogeneities lead to the acceleration of the cosmic volume 
expansion. In §3, in order to obtain quantitative information about the acceleration 
of the cosmic volume expansion, we study the stuffed Swiss-cheese model. Section 
4 is devoted to a summary and discussion. In this paper, we employ units in which 
c = G = 1 and basically follow the conventions for the Riemann and metric tensors 
and the notation used in the textbook of Hawking and Ellis. 15 ) 

§2. One scenario, a qualitative argument 

We consider an inhomogeneous dust-filled universe which is initially expanding 
everywhere. We assume that at some stage, the dust begins contracting in some 
domains and continues expanding elsewhere (see Fig. [Q) . Such a situation is consis- 
tent with our conventional picture of the real universe: Structures (stars, galaxies, 
etc.) form within contracting domains, or there might be large-scale bulk velocity 
fields due to large mass concentrations, in an overall expanding universe. Further, 
we assume that the universe is almost periodic and thus homogeneous in an average 
sense. We adopt a dust-comoving Gaussian normal coordinate system in which the 
line element is given by 

ds 2 = -dt 2 + llj (t,x k )dx i dx j , (2-1) 

where i, j, k = 1,2,3 represent the spatial components, and jij(t,x k ) is the induced 
metric on the spacelike hypersurface t = [constant], which is orthogonal to the tra- 
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jectories of dust particles. The stress-energy tensor of dust is given by 

T» v = puW, (2-2) 
where p is the rest mass density of dust and = Sq is the 4- velocity of dust particles. 




Fig. 1. An inhomogeneous universe with expanding and contracting regions. 

We consider a compact spatial domain D in each t = [constant] hypersurface 
by assuming that D is dust-comoving. By definition, D contains a fixed rest mass 
of dust. We assume that the domain D is so large that the spatial periodicity is 
recognizable in this domain. The volume Vd of D is defined by 

V D = [ ^d 3 x, (2-3) 

JD 

where 7 is the determinant of the spatial metric 7^. Following Refs. 13) and 16), we 
define the spatially averaged scale factor ao{t) of the domain by Sap /an = Vd/Vd, 
where the dot represents the derivative with respect to t. This definition is equivalent 
to the relation 

a D (t) oc v£ /3 (t). (2-4) 

Because we assume that the expanding and contracting domains coexist within 
the domain D, we rewrite the volume Vd in the form 

V D = V c + V c , (2-5) 

where V c and V c are the comoving volumes of the union of the expanding domains 
and that of the contracting domains, respectively. The time derivative of the volume 
is given by 

V D = V C + V C = \V C \-\V C \. (2-6) 

We see from the above equation that even though Vd is initially positive, once V c 
becomes dominant, Vd can be small or even negative. In other words, due to the ap- 
pearance of contracting domains in D, the cosmic volume expansion of D can greatly 
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slow down, or the volume Vo of D can decrease. However, the contracting domains 
will collapse, and then their contributions to the volume Vd will necessarily become 
negligible; in realistic situations, almost stationary structures, like stars, black holes 
and galaxies, are formed. This means that V c will be dominant in Vd, and thus the 
volume Vd begins to increase again. Here it should be noted that acceleration of 
the volume expansion is realized near the end of the slowdown or contraction stage. 
During this period, the second-order time derivative of the effective scale factor ao(t) 
becomes positive. This is the mechanism of the acceleration of the cosmic volume 
expansion proposed by Nambu and Tanimoto. 14 ) 



§3. Stuffed Swiss-Cheese Model 



In order to quantitatively study the effects of inhomogeneities in the universe, 
we construct a specific model. First, we consider the Einstein-de Sitter universe 
(EdS) and remove disconnected spherical domains from it, so as to guarantee spatial 
periodicity. Next, each removed region is stuffed with spherically symmetric but 
inhomogeneous dust whose Misner-Sharp mass 17 ) is the same as that of the original 
EdS regions. The dynamics of the inhomogeneous dust ball are described by the 
so-called Lemaitre-Tolman-Bondi (LTB) solution. We call this model the stuffed 
Swiss-cheese (SSC) model. We assume that the LTB regions are initially expanding 
but eventually begin to contract and form singularities and black holes. 

The line element of the LTB region is given by 

ds 2 = -dt 2 + ^—j^L- dx 2 + Y 2 (t, X )(de 2 + sin 2 Odtp 2 ), (3-1) 

where the prime denotes differentiation with respect to the radial coordinate x, which 
is assumed to be non-negative. This coordinate system is a dust-comoving Gaussian 
normal coordinate system, as in the previous section. The Einstein equations lead 
to the equations for the areal radius Y(t, x) and the rest mass density p(t, x) of the 
dust, 

Y' = -*(X)X 2 + (3-2) 
- - M '<*> (3-3) 



4ttY'Y 2 ' 



where k(x) an d M(x) are arbitrary functions of the radial coordinate. We assume 
that p is non-negative and Y is monotonic with respect to x, i.e., Y' > in the 
regular region. This assumption leads to M' > 0, and thus we can set 

Mix) = ^V, (34) 

where pq is a non- negative arbitrary constant. Our treatment does not lose gener- 
ality with the above choice of M{x)- Equations (|3-ip - (|3-3p are invariant under the 
rescaling of the radial coordinate x> 

X X = X(x)- (3-5) 
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By virture of this property, the above form of M{x) is entirely general. 
The solution of Eq. (|3 2|) is given by 



Y 



t - ti(x) 



47rpo 
3fc(x) 

471-10 

3A:(x) 



cos 
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sin 
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(3-6) 
(3-7) 



where ii(x) is an arbitrary function. We can use this form of the solution for any sign 
of k(x)- Note that ii(x) is the time at which a shell focusing singularity appears, 
where 'shell focusing singularity' means Y = for % > and Y' = at x = 0- 
In this paper, we consider the region satisfying t > t\, and hence the time t = t\ 
corresponds to the Big Bang singularity. Hereafter, we assume a simultaneous Big 
Bang, i.e., t\ = 0. Using the terminology of cosmological perturbation theory, there 
are only growing modes near the Big Bang singularity in this scenario. 

We denote the boundary between the LTB region and the EdS region by x = Xb- 
Then, we divide the LTB region into four regions, [0, xi), [xi>X2)> [X2,X3) an( 4 
[X3, Xb)- We use the following spatial profile of the curvature function: 



(x 2 - xi) 



Hx) 



r k 
kg 



2\2 
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+ X1 + X2 



kp 
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[xi + xi) 




for < x < Xi, 
for xi < X < X2, 

for X2 < X < X3, 

\ 2 

1 > for X3 < X < Xb, 



(3-8) 



where ko is constant. In order to guarantee the relation 1 — kx 2 > 0, the following 
inequality should hold: 



k 



(x? + xi)<i- 



(3-9) 



Since we are interested in the case that the LTB regions contract, we assume ko > 0. 

The parameter k is closely related to the value of the comoving volume in the 
LTB region Vltb ; which is bounded below as 
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(3-10) 
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where, for notational simplicity, we denote Y\ x=Xi by Y^. Note that the last term 
in the above equation, is equal to the volume of the original EdS region, 

since the areal radius Y is continuous across the boundary between the LTB and 
EdS regions. Also, the first term in the last inequaltiy gives a lower bound on the 
volume increase SV due to the inhomogeneous geometry: 



S V> ^l^^ (Yi-Yi). (341) 

In the limit k — ► 1, the volume increase becomes infinite. Thus if k is almost equal 
to unity, the volume of the LTB region can be dominant in the SSC universe, even 
if Y3 becomes small. Further, note that this result is not related to how the LTB 
region is divided into the four regions [0,xi), [xi,X2), LX2,X3) and [x3,Xb), if K is 
sufficiently close to unity. 

The physical meaning of the volume increase is the decrease of the binding energy 
in the LTB region. The binding energy -Ebind(x) within the spherical domain of the 
comoving radius \ is defined by the difference between the gravitational mass M(x) 
and the rest mass of the dust, i.e., we have 



E hi nd(x) ■= M(x) - Att / p d X , (3-12) 

Jo v 1 — kx 



where the second term on the right-hand side corresponds to the rest mass of the 
dust included within the spherical domain of the comoving radius x- in the- case of 
vanishing ko, the LTB region is identical to the original EdS region, and -E'bind(x) 
vanishes. By using Eq. (|3-3|) and the same manipulation as in Eq. (|3-1U|) . we obtain 
an upper bound on the binding energy -Ebind(Xb) of the LTB region as 

£bind(Xb) < [M( X3 ) - M( X2 )] . (3-13) 

V 1 — K 

The limit k — * 1 leads to a negatively infinite binding energy. In other words, if k 
is positive, the rest mass of dust in the LTB region is larger than that included in 
the original EdS region. An arbitrarily large rest mass can be put into the removed 
domain if k is set sufficiently close to unity. In this sense, the LTB region with k — 1 
is a very high density region compared with the original EdS region. 

The singularity formed at the origin, x = 0> can be a null naked singularity, 
whereas that formed in the domain x > is spacelike. 18 )~ 20 ) In the case of the present 
model, since the innermost region [0,xi) is a Friedmann universe with a positive 
spatial curvature, the singularity formed at x = is spacelike and thus not naked. We 
see from Eq. (|3-7p that the singularity formation time is t = 87r 2 /)o/(3A; 3//2 ). Because 
k is monotonically decreasing with respect to x, the singularity formation time is 
monotonically increasing with respect to x> an d thus no shell crossing singularity is 
realized in the present model. Therefore, singularities formed through gravitational 
collapse in the LTB regions are necessarily spacelike in the present SSC model. There 
is no causal influence of spacelike singularities on the regular domains. This property 
allows us to ignore the contribution of singularities in the calculation of the volume 
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Vd of the domain D, and thus we can follow the evolution of the cosmic volume 
expansion even after singularity formation. Further, it should be noted that the 
present choice of k(x) guarantees that the LTB region smoothly connects to the EdS 
universe. The line element of the EdS region is given by the form 

ds 2 = -dt 2 + a 2 e (t) [d X 2 + X 2 {dO 2 + sin 2 6d<p 2 )] , (3-14) 

where 

a e (t) = (6TT Po t 2 ) 1/3 . (3-15) 

We assume that the universe has the periodicity of the comoving interval, 21, 
in the original Einstein-de Sitter universe. With this assumption, a cubic region of 
comoving edge length 21 corresponds to the compact domain D considered in the 
previous section. We define the scale factor a(t) of this inhomogeneous universe as 

ao{t) = -|-. (3-16) 

3.1. One-scale model 

We first consider a model with one LTB region with Xb = I i n the domain D 
(see Fig. Ell . 




21 



Fig. 2. The compact domain D of the one-scale SSC model. 



The volume Vrj of the domain D is given by 



a c (t)i Y 2 dY 



( 4vr\ o o f^W 



V 1 - k (x)x 2 



(3-17) 



The integral on the right-hand side of the above equation corresponds to the comov- 
ing volume of the LTB region. It should be noted that x i n this integral is a function 
of the areal radius Y at time t through Eqs. (|3-6p and (|3- 7|) . Thus, this integration 
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Fig. 3. The left panel plots the evolution of the averaged scale factor od of one scale SSC model 
for various values of k. The right panel plots the acceleration of the averaged scale factor ar>. 



covers the only region which has not yet collapsed to the singularity, i.e., the domain 
satisfying Y > at time t. 

In Fig. El we plot several examples of the evolution of the scale factor ao{t). 
We set H l = y/B/2,H xi = V6/8,H oX 2 = >/6/4 and H x 3 = 3\/6/8 with H = 
^J&ttpq/'S. We chose these values so that the areal radius of the boundary of the 
LTB region is equal to the horizon radius of the EdS universe at t = I. These curves 
represent the scale factors as functions of t for various values of k defined by Eq. (|3-9p . 
As mentioned in the preceding section, the slowdown of the volume expansion or the 
contraction of the domain D actually does occur, due to the contraction of the LTB 
region. After the LTB region collapses to a singularity covered by a horizon and 
forms a black hole, the volume expansion of the domain D, which has been slowed 
by the contracting LTB region, accelerates until the speed of the volume expansion 
becomes the same as that of the EdS universe. Thus, at a time near the end of the 
slowdown or contraction stage, an acceleration period appears. In the cases depicted 
in Fig. [3l acceleration of the volume expansion occurs in four cases (k > 0.999453), 
and in these four cases, volume contraction before the acceleration period occurs in 
the two cases with k > 0.999922, whereas the other two cases with smaller k do 
not have a contraction period, although acceleration does occur. In the case of the 
smallest k (= 0.984375) depicted in Fig. [3j no acceleration period appears. 

The necessary condition of the appearance of an acceleration period is that the 
volume of the LTB region becomes dominant in Vd- As mentioned above, if the 
parameter k is almost equal to unity, the contribution of X in t ne range \2 < X < X'Z 
is dominant in the volume of the LTB region and, further, makes the volume of the 
LTB region dominant in Vo for a long time before the domain x < X3 completely 
collapses. Thus, it is conjectured that acceleration of the cosmic volume expansion 
is realized if and only if n is almost equal to unity, and we can see from Fig. [3] that 
this is in fact the case in the present model. 

The size of the LTB region relative to the cosmological horizon is very important. 
If there are many contracting LTB regions within the cosmological horizon of the EdS 
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universe in the acceleration period, the SSC model might mimic our real universe. 
By contrast, if the size of the LTB region is comparable to the cosmological horizon 
in the acceleration period, the behavior of this SSC model might be very different 
from that of our universe. To clarify the situation, we investigate which of these 
two cases corresponds to this SSC model. It should again be noted that when the 
volume of the domain \ < X3 is dominant in the volume Vr> and is contracting, the 
volume expansion of the domain D greatly slows down, or the domain D contracts. 
The acceleration of the volume expansion is realized when the volume of the domain 
X < X3 becomes so small that its contribution to Vq is negligible. In other words, 
the acceleration period appears around the time at which the dust in the domain 
X < X3 collapses to a singularity, and this time is given by 



Sir 2 poxl 

' 3 K 3 /2 ' 



At t = t ac , the ratio of the areal radius of the boundary of a LTB region a e (t)l to 
the Hubble horizon radius of the EdS region, := a e /a e = 3i/2, becomes 



Here it should be noted that k has to be almost equal to unity so that the acceleration 
of the cosmic volume expansion is realized. Therefore, the size of the LTB region is 
comparable to the horizon radius of the EdS region in the acceleration period. This 
means that non-linear general relativistic effects are important for the realization of 
the acceleration period. 

3.2. Two-scale model 

As the second model, we consider the situation that there are two kinds of LTB 
regions, one with radius Xb = I and the other with radius Xb = (v3 — 1)/- We assume 
that the large LTB region is put in the cubic domain D in the same manner as in 
the first model, whereas an eighth of the small LTB region is put on each vertex of 
D (see Fig. i]). 




(3-19) 
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In this model, the volume of the compact domain D is given by 



V D 



+ Air 



47T 

a c (t)l 







1+ h/3-l 
i - ^(xl)Xl 



a e (i)Z 3 

•(V^3-l)oe(t)J 



+ 47T 



yW 



y/i-Hxs)xl 



(3-20) 



where xl = Xl(^) an d Xs = Xs(X) are the radial coordinates as functions of the 
areal radius Y at time t in the large and small LTB regions, through Eqs. f|3-6j) and 
(|3-7|) . respectively. We present one example of the evolution of the scale factor ao(t) 
in Fig. [5j As expected, there are two acceleration periods. The first is realized when 
the small LTB regions collapse to the singularities, and the second appears when the 
large LTB region collapses. Thus, the coexistence of various scales of LTB regions 
can lead to more than one acceleration period. 
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Fig. 5. The left panel displays the evolution of the averaged scale factor for the two-scale SSC 
model. We set k — 0.999688 in the large LTB region and k = 0.999995 in the small LTB region. 
Except for these values of k, we used the same parameter values here as in the one-scale SSC 
model. 

§4. Summary and discussion 

Using the stuffed Swiss-cheese model, we investigated how the cosmic volume 
expansion is accelerated by inhomogeneities in a dust-filled universe. We found that 
the existence of contracting regions is a necessary condition to realize the accelera- 
tion of the cosmic volume expansion. Near the time at which the central contracting 
region collapses to form a singularity, the acceleration of the averaged scale fac- 
tor becomes positive, and the acceleration phase of the spatially averaged universe 
appears. Such a situation might be the same as that in the real universe, with struc- 
tures (e.g. galaxies, clusters of galaxies) forming within the contracting regions, or 
there might be large-scale bulk velocity fields due to large mass concentrations, in 
our expanding universe. However, acceleration of the cosmic volume expansion is 
realized only when the size of each inhomogeneity is comparable to the cosmological 
horizon scale in the case of the stuffed Swiss-cheese universe model. We see from 
this result that this model is very different from the universe observed today in the 
period of accelerating cosmic volume expansion. However, for the model with various 
scales of inhomogeneities, the temporal variation of the cosmic volume expansion is 
nontrivial from a theoretical point of view, and a complete analysis is left as a future 
work. The analysis of the backreaction effect for the dust-filled universe based on 
a second-order cosmological perturbation reveals no signs of acceleration of the vol- 
ume expansion. 10 )' 11 )' 13 ^ Thus, we conclude that the appearance of the acceleration 
phase in our model is due to a highly non-linear effect expected to appear beyond 
third order in the perturbation analysis. 

Here we note that the issue of cosmic acceleration first arose from observational 
data of the distance-redshift relation and that we have not directly observed the 
acceleration of the 'volume' expansion. If the universe is homogeneous and isotropic, 
then we can conclude from the observational data that the cosmic volume expansion 
is indeed accelerating. By contrast, in the inhomogeneous universe model considered 
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in this paper, the cosmic expansion of the comoving volume does accelerate, but 
it is not yet clear whether the distance-redshift relation is similar to that of the 
homogeneous and isotropic universe with accelerating volume expansion. There are 
several studies of an inhomogeneous but isotropic dust-filled universe whose distance- 
redshift relation is consistent with the observational data. 1 ) -4 ) In these models, 
there is no contracting region, and therefore the cosmic volume expansion does not 
accelerate. Conversely there is a possibility that the acceleration of the cosmic 
volume expansion due to inhomogeneities does not lead a distance-redshift relation 
similar to that obtained in the observation of Type la supernavae. Investigating this 
point is also a future work. 
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